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HYDROMAGNETIC WAVE PROPAGATION 
IN A CONSTANT DIPOLE MAGNETIC FIELD 
By John Carstoiu 


SUMMARY 


The study begins with sane formal apparatus of the theory (non-linear). 
The analogue of Clebsch' s transformation of the hydrodynamical equations 
allows a useful representation for the magnetic field; this representation 
leads in turn to some canonical equations for the motion of mag netic field 
of great theoretical interest. 

A theorem of decomposition of the magnetic field similar to the 
Cauchy-Stokes decomposition theorem in hydrodynamics is also here pre- 
sented. The rest of the report deals with hydrcmagnetic wave propagation 
in (i) a constant and uniform magnetic field; (ii) in a constant dipole. 

The effect of compressibility is especially investigated in the case 
of a constant and uniform magnetic field. 

The disturbance is specified in terms of vorticity and current 
density. It appears that the compressibility of a medium acts as a wave 
filter discriminating between components of vorticity (and current 
density) and passing only those directed along the (undisturbed) magnetic 
field. 

The case of a dipole magnetic field presents a singular importance, 
in view of its applications to geophysical phenomena, and is discussed in 
some great detail. Dungey in his remarkable report of 1954 (The 
Pennsylvania State University, Ionisphere Research Laboratory, Scientific 
Report No. 69 ) has already discussed the electrodynamic behavior of the 
Outer Atmosphere in the presence of a constant -dipole magnetic field. In 
this study, however, the problem is approached from a different point of 
view; the magnetohydrodynamic behavior of the fluid is discussed in terms 
of vorticity and current density. The equations obtained are complicated, 
however, and solutions are discussed only at large distances from the 
center of the dipole. 

The study of hydrcmagnetic wave propagation in a dipole is preceded 
by a chapter where the geometry of lines of force is presented. 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-1689 

HYDROMAGNETIC WAVE PROPAGATION 
IN A CONSTANT DIPOLE MAGNETIC FIELD 
By John Carstoiu 

Sylvania Electric Products, Inc. 
Waltham, Massachusetts 


INTRODUCTION 


The theory of magnetohydrodynamic waves is the creation of Alfven 
and Walen. In his celebrated paper of 1942 [(l)J Alfven has shown that 
waves can travel along magnetic lines of force in a conducting material; 
see also Reference 2. Walen, [(24)] in 1944, has discussed in some de- 
tail these waves and has given the magnetohydrodynamic equations starting 
with the principle of conservation of energy. 

To these pioneer efforts other savants (Spitzer [( 19)1 > Cowling C(7> 
8)1, Grad [(15)] , Lighthill [(17)1, and MacDonald [(l8)]) added their 
researches, resulting in the elegant theory of magnetohydrodynamic waves 
as presented, for instance, in Lighthill ’s great memoir of 1961. Ity 
earlier work consists mainly in analysis for the case of a compressible 
fluid [(3,5)] • 

All this work presupposes a uniform magnetic field. Non-uniformity 
of the magnetic field affects the theory both because the wave velocity 
varies in magnitude and direction and because new forces are introduced. 
The case of a dipole magnetic field is of particular importance^ in view 
of its applications to geophysical phenomena. This revivifies Stormer's 
work [( 22 )j and other aurora theories. It is of great importance indeed 
to know up to what an extent Stormer’s theory may be improved by using 
the hydrcmagnetic approach, that is to say, assuming a fluid mechanics 
continuum approximation. 

hydrcmagnetic wave propagation in a dipole magnetic field has 
already been discussed in seme detail by Dungey [(9)] • approach is, 

however, more general and sets up the magnetohydrodynamic equations in 
the Outer Atmosphere at large distances from the earth. 


In a subject which is developing so rapidly, this problem cannot be 
discussed with any approach to finality, but this did not seem to be a 
reason against writing these pages. 


In Chapter I we rapidly review the basic equations of magnetohydro- 
dynamics. . The analogue of Clebsch's transformation of the hydrodynami- 
cal equations allows a useful representation for the magnetic field; 
this representation leads in turn to sctne canonical equations for the 
motion of magnetic field of great theoretical interest. The Chapter is 
concluded with a theorem of decomposition for the rate of change of the 
magnetic field similar to the Cauchy-Stokes decomposition theorem in 
hydrodynamics. Chapter II presents the theory of magnetohydrodynamic 
waves in the presence of a uniform magnetic field. In Chapter III we 
discuss in sane detail the geometry of dipole magnetic lines of force. 

Chapter IV deals with magnetohydrodynamic waves in a constant dipole mag- 
netic field. 6 


The material in this report summarizes the results of researches 
undertaken by the writer on the subject of Hydromagnetic Wave Propagation 
in an electrically conducting fluid of infinite extent embedded irTa con- 
stant dipole magnetic field. This work has been supported by a contract 
MAC , National Aeronautics and Space Administration (Contract No. 

NASr- 18 ) The little time allowed to us to investigate the subject and 
the complexity of the problem have permitted only an exploratory effort 
The phenomena, as indicated by mathematical analysis, are no doubt very 
complicated and considerable more effort is required in this direction 
It is hoped, however, that this pioneer effort may serve as an introduc- 
tion to subsequent detailed investigations. 

t + Ih 0 author takes this occasion to express his gratitude to Dr. Robert 
Jastrow for his interest in and support of this work. 
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CHAPTER I 


BASIC CONCEPTS - THE EFFECT OF A MAGNETIC FIELD 


1 . The Basic Equations of Magnetohydrody namic s 

When the displacement currents may be neglected, Maxwell’s equations 

are 


curl H = tatJ y 
c>H 

curl E = -P e ^ , 


( 1 ) 

( 2 ) 


div H = 0 , 


(3) 


where the electranagnetic variables are measured in e l ect J^ et i C .^ J 
E and H are the intensities of the electric and magnetic fields, J is th 
current density, and p e is the magnetic permeability. To complete the 
equations for the field, we need an equation for the current den y. 

Consider an electrically conducting fluid which has a conductivity a 
and executes motions described by the velocity v. The elec nc le i 
will experience is E + p e v x H, thus 


J = cr(E + 


P e v x H) 


w 


The equations (l) - (*0 incorporate the effect of fluid motions on 
the electranagnetic field. The inverse effect of the field on the motions 
results from the ponderonotive which the fluid elements ^erience by 
virtue of their carrying currents across magnetic lines of for . 
the Lorentz force given by 


= p J x H = curl H x H . 


(5) 


Including this force among other forces acting on the fluid, we have 
the equation of motion 


dv 

p — = div P + pX + p J x H , 
at — e 


( 6 ) 
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where p is the density ]? is the total stress tensor and X represents the 
external forces of non-elect ronagnetic origin. 


In tensor notation^ this equation can he written 


P 



+ v 



H dH. dp. . 

5 l H j S 3 T = ^ + px i 
j j 



where explicitly 


ij 


= -p 5 


ij 


+ 2pe. . 
ij 


2 . 

- :=■ p6. .e. . 
3 ij kk 


( 7 ) 


(8) 


and. where p is the isotropic pressure, p is the coefficient of* viscosity , 
and ej_j is the rate of deformation given by 



For an incompressible fluid in which p is constant and the forces X 
derive from a potential - fl, the equation of motion (7) simplifies to ~ 


8 v. by. p H. 8 h -s / p h 2 \ 

■ST + v j - 10 SJ = - SH* 1 + P + TS-] + >'v 2 v i , ( 10 ) 

where v = p/p denotes the kinematic viscosity . 

In the general case, the equation of motion (7) has to be supple- 
mented with the equation of continuity 


8p 8 , N 

3t + 577 (pv j } " 0 > (n) 

J 

and the heat equation. We shall not write down the heat equation assum- 
ing in this study that our variables do not depend of temperature 


k 



2. The Equation of Motion for the Magnetic Field 

We shall now obtain an equation of motion for the magnetic field. 

In view of further developments, it is convenient to introduce here the 
vector potential A and the electrostatic potential 0, writing in the usual 
way 


H = curl A , 


div A = 0 , 


3a 

E = - - grad 0 . 

We then obtain according to equation (4) 


J = a 


3 A 

5t 


- grad 0 + p g v x curl A 


Substitution into equation (l) gives 


( 12 ) 

(13) 

(14) 


(15) 


3a 

-r— = v x cirri A - grad 0 - v curl curl A , 
ot — — e 


(l 6) 


where v e = (W e a) -1 will be designated as the magnetic viscosity (see 
Elsasser E(l3)3, page 21 ) . It may be noted that v e like v is of dimen- 
sions cm2 sec -1 . 

Taking the curl of terms of equation (l6) and assuming v g a con- 
stant, we obtain (in Cartesian coordinates) 


3h 
5t = 


curl(v x H) + v 


(17) 


which is the equation of motion governing magnetic field. Equation (17) 
is general; it is not restricted either to incompressible fluids or to 
inviscid fluids. 

The case when the electrical conductivity of the medium may be con- 
sidered as infinite is a particular interest in cosmic electrodynamics. 
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The magnetic viscosity is then zero, and equations (l6) and (17) reduce 
respectively to 

8a 

^ = v x curl A - grad 0 , (18) 

8 h 

= curl (l x H) > (19) 


Equations (18) and (19), especially the latter, have been the object 
of considerable research in the literature. We may notice at once its 
full analogy with the Helmholtz equation for the vorticity. This immedi- 
ately permits to apply mutatis mutandi s the classical and elegant results 
of the theory of vorticity to the magnetic field. In particular, it 
follows that the lines of force move with the fluid . For further details 
of this analogy we refer to Goldstein's Lectures on Fluid Mechanics I”(l 4 ) 
page 76. L - 

3 * The Elsasser-Carstoiu Theorem 


Equation (19) may be put in a form which is reminiscent of the 
Cauchy-Stokes decomposition of an arbitrary instantaneous continuous 
motion of a fluid (see for instance Truesdell f(23)] , page 65). Equiva- 
lent to the basic equation (19) is 


_d_ 

dt 



H. dv. 

_J i 

P dx. * 
J 


( 20 ) 


where the equation of continuity (ll) has been used. Equation (20) can 
be rewritten 


_d_ 

dt 




(21) 


where, besides the rate 


of deformation e-j_j. 


the 


vorticity 0) = tujj 


to 


ij 



(22) 


appear under its tensor components. We can write 
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(23) 


d 

dt 



H. 

= — J- a> 

P U 


1 8g 

+ 2?ShT ' 

1 


or in vector notation 



— i 

(i) x — + ~ — Grad.- G , 
- P 2p H 


(24) 


where we set 


e , .H.H. 
ij i J 


(25) 


and the gradient is taken with respect to H^. 

Equation (24) shows that the rate of change of the magnetic field 
niny he conceived as made up of two parts. The first part expresses a 
rotation of the field with the fluid particle; the second part shows that 
the terminus of H is moving in the direction of the normal to that 
quadric of the system : 

e. .X.X. = const. (26) 

ij i J 


on which its terminus lies. 

In this form the theorem has been stated by the writer [(4)] . An 
integral formulation closely related to this has been earlier given by 
Elsas ser ^(11,12)]] . 

4. The Analogue of the Clebsch Transformation 

A matter of interest in cosmic electrodynamics is the analogue of the 
Clebsch transformation of the hydrodynamical equations (see Lamb Ql6)J, 
page 248). Putting 

A = grad P + S grad i|r (27) 


one has 

H = curl(!E grad \|r) = grad $ x grad \|r 


(28) 


7 


The representation (28) is identical to that given by Sweet [(2l)] 
(see also Dungey [(lO)j , page 31 ) with the exception of a factor F, func- 
tion of 0 and \|r only, which appears on the right side of ( 28 ) in Sweet's 
representation. However, it can be shown that Sweet's and our repre- 
sentation are equivalent (see Lamb, loc. cit. ). 

The immediate consequences of formula (28) namely 


H • grad S = H • grad i|r = 0 


(29) 


show that the magnetic field is tangent to the surfaces SB = const, and 
= const., which we shall call surfaces of force , and which correspond to 
vortex surfaces in hydrodynamics. It is evident that their intersections 
are the lines of force. 

5- Hamiltonian Form for the Equation of Motion for the Magnetic Field 

Let us now cane back to equation (l 8 ) and substitute thereto the 
value of A given by (27). We have 


(grad P + 5 grad t|f) = v x (grad $ x grad \|r) - grad 0 

= (v • grad \|r)grad 5 - (v . grad B)grad 
- grad 0 , 


♦ 

(30) 


which can he written 


d® 

dt 


grad ^ 


df 

dt 


grad & = -grad pC 


where 


jc = 


8 p 

St 


+ s 



(31) 


(32) 


Scalar multiplication of terms of equation (3l) by grad 3 ) x grad \|r 

gives 

grad/^ • (grad $ x grad i|r) = 0 , ( 33 ) 

that is the Jacobian 


8 



This shows that$> is of the form^(!E^,t) . Hence 


graul H = ^ grad t + ^ grad fi . 

Comparison of equations (3l) and (35) gives at once the Hamiltonian system 

d® cKf = ^_ ( 36 ) 

dt 5^ * dt 5® 

Equations (36) are analogous to Stuart’s equations in hydrodynamics (see 
Lamb, loc . cit.) and were derived by this writer in a recent paper l^J * 


CHAPTER II 


MAGNETOHYDRODYNAMIC WAVES IN A CONSTANT AND UNIFORM MAGNETIC FIELD 


6. The Case of an Incompressible Fluid 

We begin with a discussion due to Walen [(2^)] (see also Cowling [(7, 
8)J . Consider an infinite mass of uniform fluid, at rest, embedded in a 
constant and uniform magnetic field H 0 . We assume this fluid to be an 
inviscid, incompressible and perfectly conducting material (c— >co). 
Suppose that as a result of a perturbation, a velocity field v is pro- 
duced in a certain region, and that the magnetic field becomes Hq + h. 

The equations giving the variations in v and h are 


3v 

3 o St 


p 


= - grad p + p^g + curl h x (H^ + h) 


(37) 


Sh 

St = curl 


V X 


(Ho + h) 


(38) 


where we have included, the gravitational potential 


g = - grad ft 


(39) 


P Q is the uniform density of our fluid, and the term (v -V) v has been 
emitted, ” 


Now, since H = constant and 
7 — o 

div v = 0 , 


(to) 


div h = 0 , 


equations ( 37 ) and ( 38 ) simplify to 


8v 

St 


grad 


P + 


P H 
e— o 


4it 


p o n| 


+ r^- (H 

4TC ~ O 


V)h 


Sh 

St = (So * v) - , 


(4l) 


(42) 

(43) 
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by neglecting squares and products of the small quantities h,v. 
the divergence of equation (42); we have 

/ pH • h 

div grad | p + ^ + P Q 

In Cartesian coordinates equation (44) becomes 



Take 


(44) 


J p e»o ■ * 

v Ip + 


bn 


4“ P ^ — 0 y 

o 1 


( 45 ) 


since p + (p • h) /kit + has no singularities and is bounded, 


p H • h 


P + 


e— o 


t + n H = constant 

bn o 


( 46 ) 


Hence, equation (42) bee ernes 


4jtp yr — p (ll • V)h 
K o ot e — o — 


( 47 ) 


For simplicity take Oz parallel to H^. Then equations (47) and (43) 
become 


bv bh 

4jtp o si = P e H o ’ 


Sh Sv 

5t _ H o 3^ 


Hence, by cross differentiation 



( 48 ) 

( 49 ) 


(50) 


ll 


d 2 h 2 d 2 h 

a? " A ° a? ' 


where 



( 51 ) 


(52) 


is the Alfven’s phase velocity, named so in honor of its discoverer. Thus 
the disturbance can be expressed as the resultant of two sets of waves 
traveling with velocities +A 0 in the z-direction, i.e., along the lines of 
force of the undisturbed field m — - -- - 

(m-h. ) waves . 


These waves are called magnetohydrodynamic 


After the two waves have separated we have in either of the waves 


dv dv 

= ±A o ' 


(53) 


the sign depending on the direction of propagation 
Comparison of equations (49) and ( 53 ) gives 


of the wave considered. 


H v 




(54) 


a g ° in &, farther > we no- te, that in considering the propagation 

described by Alfven, the velocity v and the magnetic field h can be replaced 
by the vorticity ca and the current density j = (l/ 4 jr)curl h respectively 
for one has similar equations for these quantities, namely - 



(55) 

(56) 
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together with the relation 


i = 


H 

+ 

- 2rtA 

o 


CD 



( 57 ) 


7. Compressible Fluid - Vorticity and Current Density Propagation 

In taking this point of view as a point of departure, we shall show 
that in the case of a compressible medium the components of cd and in the 
direction of the field only are propagated in Alfven's manner . Thus , . 
surprisingly enough, the compressibility of a medium acts as a wave filter 
discriminating between components of vorticity and current density and 
passing only those directed along the (undisturbed) magnetic field. The 
proof goes like this. When the compressibility is taken into account, the 
linearized system replacing equations (42) and (43) is 


Sv d H 8h 

p o si = - grad * + "fe 2 5 

Sh 8v 

St = H o S^ ' H o div i » 

where we set 

d H . h 

A 0“O — n 

0 = P + 4 — + P G fi • 

Equations (40) and (4l) are replaced by 


(58) 


(59) 


(60) 


Sp 

St 


+ p o 


div v = 0 , 


(61) 


div h - 0 , 


(62) 


where p is the perturbation in density. We shall assume that 


P 



(63) 


where a is the ordinary sound speed in the absence of a magnetic field, 
o 
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Taking the curl of terms of equations (58) and (59) we obtain 


ocn 

2p o St = P e H o 5z 


(64) 


Sj_ S(i3 

2+71 5 t = 2H o Sz + x grad div 


(69) 


which imply important consequences, as will be shown. 

(a) Propagation of z-Ccmponents of Vorticity and Current Density 
Equations ( 54 ) and (65) when projected on the C> z axis give 


So3 


Sj r 


2p o "ST = P e H o “Sz" 


(66) 


2jt 


FT = H o 


Sa> 

z 

dz 


(67) 


Hence 





St 2 




and 


J 


z 


H 

+ _ 2 _ 

“ 2 rtA 

o 


O) 

z 


( 68 ) 

( 69 ) 


( 70 ) 


Thus, the components of 03 and j along the lines of force (longitudinal 

components) are propagated in the opposite directions of the undisturbed 

field with velocities +A . 

— o 


Ik 



The coupling relationship ( "JO ) between longitudinal components shows 
that (i) it does not depend on the magnitude of the magnetic field pres- 
ent; (ii) the vanishing of either component involves the vanishing of the 
other; this occurs when either quantity is zero initially. 

(b) Equations for the Transverse Components 

It may be noted that although uw and j z are propagated one -dimensionally, 
along the mag netic lines of force, no other component of vorticity and cur- 
rent density is; the x- and y-ccmponents of vorticity and current density 
satisfy 



x 


St 2 



.2 Q 
^o d 3 p 
2p Q dtdyx3z ’ 


( 71 ) 



*3 P 

2p dtdxdz , 
o 


( 72 ) 


and 



iii 

4 * P o St 2 Sy ' 


( 73 ) 



H o a 3 P 

3t 2 Sx 


( 7*0 


showing that only for incompressible flow do they satisfy the same equa 
tions as and j z ; in a compressible flow their oscillations are coupled 
to those of density. 

8. Wave-Motion Equations for the Density and Transverse Components 
Differentiation with respect to t of equation (25) gives 




A 2 p 
o o 


H 

o 



( 75 ) 
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where the equations (58) and (6 3) have been used, and the gravitational 
potential has been emitted. Now 




A - - -£] • 


Hence, the density satisfies 



2 2 

a o v P 


knA 2 p 
o o 

H 

o 





(76) 


( 77 ) 


To obtain an equation for p alone, we eliminate j x and between equa- 
tions (73)^ (7*0; and (77)- The result is 


&L/af p 

8t 2 y8t 2 








.( 78 ) 


where a two-dimensional Laplacian appears. This equation shows plainly 
the radical departure of the "new sound-wave equation" from the ordinary 
sound equation in the absence of a magnetic field. 

We may ask now if it were possible to satisfy both equations (78) 

and 




(79) 


Then, by virtue of the former, we must have 



8 x c 



> 


(80) 


and hence equation (79) reduces to 


afp 

St 2 


2 Sjd _ 
o v 2 “ ° ’ 


(81) 


which admits a solution of the form 

1 6 



provided that 



( 83 ) 


We can now easily verity that (82) is effectively a particular solu- 
tion of equation (78) under condition (83)- Thus sound waves appear 
possible in a conducting fluid penetrated by a uniform magnetic field, 
with this great difference that they do not spread out three -dimensionally 
as in ordinary acoustics; instead, they propagate (without attenuation) 
one -dimensionally, along the magnetic lines of force. It is also inter- 
esting to note that in contrast to m.h. waves, this propagation does not 
depend on the magnitude of the magnetic field present. 

Equation (78) can be rewritten as follows 


W&fg. 

St 2 y8t 2 


(a 2 + A^V 2 


\ 


+ a 2 A 2 V p 

° ° 8z 2 


V^p = 0 . 


(84) 


Differentiation with respect to t of terms of equation (84) yields an 
equation given by Lighthill £(1?)] for the expansion A = div v. 

Let us cane back to equations ( 71 ) - (74) and improve our results. 
Elimination of p between equations ( 71 ) and (84) and then between (72) and 
(84) gives 


Jk 

(at 2 


8 2 / 2 , .2x2^ 

' (a ° ° V 


2.2 

+ a A 


-T> 


.2 


0 


-2- ■£_ - (a 2 + A 2 >7 2 ’ 

st 2 U 2 ° ° 


2,2 
+ a A 
o o 


V 


03.. 


I 03 


= 0 . 


(85) 


Similarly, elimination of p between equations (73), (74) and (84) gives 


a 2 

^ o 

S 2 fh 2 , 2 2 

p( p " ( a + A )v 

2.2 8 2 J 2 . 

+ a A — — v 


ot 2 \8t 2 0 0 i 

— V / 

oo^2 
' oz 



-l 


f- 

|at 2 

, 2 2v 2^ 

- (a + A )v 
O ° j 

I 2.2 8 2 2 

+ a A — 7 r v 

1 0 °a 7 

"N 

► 




J 


,2 

0 Sz 2 [^t 2 


Equations ( 49 ) and (50) show that the quantities 


= 0 


(86) 


/: 


L a t v 


~ - ( 

8t 2 


a 


+ A - )v ' ♦ 4*0 ^ V 2 


0) ,00 

x y 


Jyj j 

x y 


(87) 


are propagated along magnetic lines of force at Alfven velocity A_. These 
quantities are identical zero if they were zero initially; under this con- 
dition, equations (85) and (86) reduce to fourth-order equations of the 
same type as equation ( 84 ). 
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CHAPTER III 


THE GEOMETRY OF DIPOLE MAGNETIC LINES OF FORCE 


9- Preliminaries 


The equations of motion in the case of a constant magnetic dipole are 
complicated and require some detailed discussion of the geometry of the 
magnetic lines of force. Stormer [(22)] was the first to use in his equa- 
tions the arc s of a line of force, instead of the time (see his book, 
page 215). He also calculated the radius of curvature of the trajectories 
in the equatorial plane of the dipole (loc. cit., page 22l). Also, other 
authors used at least the element of length of magnetic lines of force in 
various of their calculations. However, nowhere do we find a detailed 
discussion of the geometry of these lines. In this Chapter, we propose 
to fulfill this need by a systematic account of the geometry of magnetic 
lines of force for the case under consideration. 

10. Equations of Lines of Force and Linear Element 

As well known, a dipole magnetic field has components 


H 

x 



H = -M 2X2 

y ^ 



(88) 


2 2 2 2 

where M is the magnetic moment of the dipole , r = x + y + z , and the 
sign is chosen such that H z is positive in the x,y plane, which is the 
equatorial plane of the dipole; this requires that the dipole has its 
negative pole directed upward. 

The differential equations of lines of force are 


dx _ dy dz 

- 3 xz “ - 3 yz “ 2,2 0 2 

^ x + y - 2 z 

The first two equations give at once 


y = Cx , 


(89) 


(90) 


where C is a constant. Equation ( 90 ) represents a family of planes 
through Oz axis. 
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Then, by virtue of (90), the last two equations give 


2 2 2 
dz _ 2z - PCrf 

dx 3xz ’ (91) 

A 2 = 1 + C 2 . ( 92 ) 


This is a honogeneous equation of the first order, 
z = xu , 

the equation becomes 


On substitution of 

(93) 


dx 

x 


3udu 


u 


+ A 


(94) 


in which the variables are separated; the solution is 

_ 2 

x 3 = B(u 2 + A 2 ) , 

B being a constant of integration. Hence 

X 3 = B(x 2 + y 2 + z 2 ) , 

that is 
2 

= ar , 
where a = Jb . 

We now introduce the polar coordinates 

x = r cos 0 
y = r sin 0 
z = r sin X 


cos X , 
cos X , 


(95) 


(96) 


(97) 


(98) 
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where X designates the magnetic latitude. In these coordinates, 
equation (97) becomes 

2 2 3 

cos 0 cos X = a r . 

To eliminate 0, we write 

^ = C = tan 0 . 
x 

Hence 

2 ri 1 

cos = -5 > 

C + 1 


( 99 ) 


( 100 ) 


( 101 ) 


and equation (99) becomes 
2 . 

r = r cos X , 
o 


( 102 ) 


where we put r = l/[a 3 (C 2 + l)] . Equation (102) is the equation of lines 
of force in each meridian plane j it is obvious that is the value of r 
for X = 0 (in equatorial plane). 


In Cartesian coordinates, we have the following parametric equations 
of the line of force 


r cos X 
x = > 


"N 


v/c 2 + 


Cr cos X 

y = - > 

/c 2 + 1 

z = r sin X . 


J 


(103) 


To calc ula te the linear element of these lines we may use either equation 
(102) and then 
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(104) 


,2 ,2 2,2 

ds = dr + r 6.X , 

or the parametric equations ( 103 ) and then we have 
2 2 2 2 

ds = dx + dy + dz . ^ 105 ) 

The result is 

1 

o p 

ds = rdl(l + 4 tan X) . (106) 

11- The Prenet Formulas for a Line of Force 


The tangent to a line of force is defined by 
dx 3 sin X ''N 


a = 


ds 


J* + 1 Jl 


P = 


= _ 


ds 


/ 


1+4 tan X 


sin X 


C 2 + 1 Jl + 4 tan 2 A. 


3C 


7 = 


o 

dz (l - 2 tan X)cos X 


ds 


/T7 


4 tan X 




(107) 


We next calculate the quantities da/ds, dp/ds, d?/ds. 
instance, 


One has for 


da _ da dX 
ds dA, ds ’ e c ' 


After sane calculation, we obtain 
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da 

ds 


2 

3 cos X 1-4 tan X 
y c 2 7 1 r(l + 4 tan 2 X) 2 






2 1 

cL3 _ -3C cos X 1-4 tan X I 

ds ” JF71 r(1 + 4 tan2 x)2 ’ 


d7 

ds 


-9 sin \(l + 2 tan X) 

— 2 2 

r(l + 4 tan X) 


The radius of curvature p is given by 



9(1 + 2 tan^ \)^ 
r 2 (l + 4 tan 2 X) 3 


In the equatorial plane 
r 

r o 
p - 3 - 3 • 

As a verification of result (109), we may use the formula 


P = 


(1 + r' g ) g 


2 2 
r + 2r' 


rr 


One has 


(108) 


(109) 


(no) 


( 111 ) 
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r 


2 

= r cos 
o 


\ , 




r 1 = -2r cos X sin X , 
o 9 

r" = -2r (cos^ X - sin 2 X) . 
o ' 


(112) 


Hence 


3 

r^cos^ 1 X + 4 cos 2 X sin 2 X) 2 

2 I” 2 2 2 o "t 

cos X|_cos X + 8 sin X + 2(cos X - sin X)J 

3 

2 2 2 
r cos X(l + 4 tan X) 
o 

3(l+2 tan 2 X) 


(113) 


which checks our earlier result (109). 

principal normal, which in our case reduces to the normal of 
lines of force as these are plane curves , is determined by 




a i = p dF 


cos X 1-2 tan X 


Jc 2 + 1 Jl + k tan 2 X 


o CcosX 1-2 tan X 

1 ds R / — 

yC + 1 yj 1 + 4 tan X 




dr 

r i ■ p s 


I 


3 sin X 


1 + if tan X 


J 


(H4) 


The hinormal of lines of force is, of course, the unit normal to the 
planes y = Cx; therefore, its direction causes are 


24 



a 2 = 


- — 7 , 3 2 = - 7= 

7c 2 + i 7c 2 + i 


> 7 0 = 0 • 


(115) 


Therefore, direct application of the Frenet formulas has to give the 
above results. This will verify our previous results. One has 


a o = = , ; L 

7c 2 + i 1 + 4 tan x 


[ * P o 2 \2*1 

9 sin X + cos X(l -2 tan X) j 




C ^ — (1 + 5 tan X + 4 tan X) 


2 ,1+4 tan X 

C + 1 




2 1+4 tan X 

0 ' 


£2^ — ^ — (1 + 4 tan 2 x)(l + tan X) 


(116) 




C 2 + 1 


Also 


P 2 = 7«1 - a7 x = - 




+ 1 


r 2 = a?! - POp “ 0 » 


(117) 

(H8) 


•which values agree with those given by (115). 
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CHAPTER IV 


MAGNETOHYDRODYNAMIC WAVES IN A CONSTANT DIPOLE MAGNETIC FIELD 


12. Preliminaries 


Dungey [(9)] in his remarkable report of 1954 has discussed in some 
detail the electrodynamic behavior of the Outer Atmosphere in the pres- 
ence of a constant dipole magnetic field. We shall here approach the 
problem from a different point of view, however, concentrating our atten- 
tion, as we did in Chapter II, on the vorticity field and the current 
density. 

Consider an infinite mass of an electrically conducting fluid at 
rest embedded in a constant dipole magnetic field H. To simplify the 
discussion, take the conductivity as infinite and assume the fluid to be 
a homogeneous incompressible material. Assume that as a result of a per- 
turbation, a velocity v is produced in a certain region and that the 
magnetic field becomes K + h. The amplitude is assumed to be small 
enough for non-linear terms to be neglected. We propose to investigate 
the magnetohydrodynamic behavior of the fluid in terms of generalized 
Alfven waves. 

13- Fundamental Equations 

The relevant equations for the problem are 


Bv 


p o St = ~ grad p + v x ^’ 

(119) 

Bh 


= (H • v)v - (v - V)H , 

(120) 

div v = 0 , 

(121) 

div h = 0 , 

(122) 


the condition BH/Bt — O(constant dipole) has been used in equation (120). 
Since curl H = 0, equation (119) can be rewritten as follows 


By 

p o 5t = _srad 




(123) 


26 



Taking the curl of terms of equations (123) and (120) we obtain 


)i r -1 

2p Q = ]j^ curl [(H • V )h + (h • V )HJ , 

(124) 

. 8j_ _ _ 

4it = curl|^(H • V)v - (v • V)H 

(125) 

We have 



curl[(H • v)h + (h • v)H 

= 4jt(H • V) j + grad h 

8h 8h 

x 5J + srad h y x gj 


8h 

8h 8h 


+ grad h,x^-+grad H x x^-fgrad H^xgy 
bh 

+ grad x . (126) 


After same calculation, we obtain 

curl £( H • V)h 4* (h - v)l^ = * V)j_ ~ (j_ * V)^j > (127) 


where the condition curl H = 0 has been used. 
Hence equation (124) can be written 

2P 0 = P e [(« ■ V)j - (i • V)H] . 

On the other hand, 


(128) 


bv dv 

curl^(H • V)v - (v • V)H] = 2(H • v)u> + grad H x x ^ + grad H y x ^ 


bv 

+ grad H z x ^ 


8h 

grad V X x 5^ 


8h 8h’ 

+ grad v y x ^ + grad x ^ 


(129) 
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In this case, there is, however, no simple way to write vectorially 
equation (125) in a compact form; we have for components the following 
equations 


2 it 


■3F = [ ( ^ • v) ^Jx + ^7 e 3! + e 32 + -5# e 33 


SH 8h 

x . y 2 

e 2i + ~st e - + e - 


J 22 + e 23. 




c>H SH Sh 

x , y 

"5T e n + 


[(H • V )a> + ~xr e n + -X? e 10 + e. 


12 "5T e l 3 


'8h 5h c>H 

x , y 2 

35T e 3 i + -ST e ™ + e ' 


'32 ^x - 33 


2it 


"St" 


r 1 SH SH SH 

(H • v)uq + -5-2 e , .. + -5-^ e + -r-2 e 

L - -J z dx 12 5x 22 Sx e 23 


'Sh Sh 8h 

■37 e ll + ^ e 12 + 157 e 13 1 


(130) 


where e ± j is the rate of deformation and where again this condition 
curl H = 0 has been used. 

1^- Propagation at Large Distances 

These equations are rather complicated. We may simplify them by 
observing that the derivatives of the components of the dipole magnetic 
field are the order of r - ^ while these components themselves are of the 
order r-3. Therefore for r sufficiently large, we may neglect the term 
(j. ‘ V)H in equation (128) and similarly we may neglect all terms such a 1 
(bE x /by)ej 1 , etc., in equation (130). Hence, for r sufficiently large, 
equations (128) and (130) reduce to 


So> 

2p o 57 ■ t*- lH 


/ S A °i aj\ 

ey X 57 + H y 57 + H z 3i :) ’ 


(131) 


and 
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( 132 ) 


deo den den 

2n — = H t — + H t — + H t— • 

St x dx y oy z dz 

In the second place, we have a long a line of force 


H x = H S> *V = H S’ H z = H i- 


(133) 


where H is the magnitude of the dipole magnetic field and ds the element 
of length of line of force. Hence along a line of force, equations (131) 
and (132) simplify to 


*5 foi dx ^ dy ^ ■ 

2p o St = ^e^Sx d^ + 5^ ds + SI dFj 


(134) 


and 




/s 


2 * St 
that is 


= H 


'CD - V^UL» T I 

-dx - d£ - dz \ 

5x ds 5y ds Sz^ ds [ ’ 


den 


den 


(135) 


den dj_ 

2p o 5i = P e H Si ' 


(136) 




dcu 


2 * St “ H Si" * 


(137) 


Observing that dH/dt 
differentiation 


= 0 (constant dipole), we obtain by cross- 


d 2 en 


dt £ 


,a fa 1 a£ *2 

“ A ^22 ds SIT ’ 


(138) 




.2 ^ A . 1 dA‘ 


_ 2 

5? = A “ i? + s ^ ^ ' 


( 139 ) 


where 
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0 pH 2 p M 2 2 , /, , 2s 

A 2 _ _e _ e cos X(l + 4 tan X) 

4:rtp Wp 5 

o o r 

2 1+4 tan 2 X 

7 6 10 " ’ (i4o) 

r cos X 


where we put 

2 ^ 

7 4jtp o * (l4l) 


We achieve the reduction of these equations by taking instead of s 
the magnetic latitude X as independent variable, We have 


cku bad j _ 

— — aX 

5s* 5\ ds 9 


(142) 



fWa if 


ckO 

dX 



and 



d 

/ah 

dX 

WJ 


dX 

ds 


( 143 ) 


(144) 


Now (see Chapter III, equation (106)) 
dX _ 1 

ds i > (145) 

r Q cos 2 X(l + 4 tan 2 X) 2 

therefore, we have 
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( 146 ) 


d 2 A. _ 2 tan X(2 tan 2 X - l) 

ds 2 r 2 cos X(l + 4 tan X) 


2 • 


On the other hand 


2 2 
dA dA dX 


ds 


dX ds * 


and 


dA 2 ^ dA 2 kx\ 2 


8 co 


ds 5s" dX yds J 5x 

67^(3 + 8 tan 2 X)tan X 

= “ n tit " I 2 7 5x 

r cos \(l + ^ tan X) 
o 

Substitution of these values in equation (138) gives 


(147) 


(148) 


S 2 cu 

p 

21+4 tan X 

cv 1 

ro 

IS 

1 

St 2 

7 “6 IcTT 

r cos X 

0 j 

Sx 2 r 2 cos^ X(l + 4 tan 2 X) 
0 


2 tan X(2 tan 2 X - l) 

^-2 3 tan X(3 + 8 tan 2 X) 

+ x— 7 —rr- r 1 ^ 


r cos X(l + 4 tan X) 
o 


(149) 


that is 


>v2 

o co 


St" 


"B ITT 

r cos X 
o 


^2 
O 0) 


Sco 


— + 7 tan X • ^ 

Sx 


(150) 


and a similar equation for j_. Supposing that 
o)(X, t) = Ae^°^ CD-, ( X) , 


( 151 ) 


we get 
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(152) 


d 2 ox 


-1 


„2 8 14 


+ 7 tan X -r— + 0 r cos X • m =0 
d aA, O —1 


d\ 


2 2 2 
where (3 = a // . 


15 • Integral Equations for Vorticity and Current Density 

Equation ( 152 ) can be transformed in an integral equation similar to 
that given by Dungey (see Reference 9 , page 33). In order to do this we 
use the identity 


^[ sec x A ( ^L sec3 X) ] = sec ^ 


-2 

d m^ 


LdX 

+ 3(l + 5 tan 2 

Hence equation ( 152 ) can be written 


dm 




(153) 


^[ SeC X A ( -l Se ° 3 X) ] = [ 3 ^ + 5 tan 2 X) 


sec X 


28 13 ,1 3 , 

X| sec X 


P r cos 
o 


(154) 


or, by putting fi = m sec 3 X , 

^[ Sec f^J = |^3(l + 5 tan 2 X)sec X - 0 2 r® cos 13 . 


(155) 


Assuming that for X = 0 , dtt/dX = 0 , i.e., dm /dX = 0, we obtain the 
following integral equation 1 


X X' 

n(\) = 3(0) + J cos V dX' J ^3(1 + 5 tan 2 X")sec X" 


3 2 r 8 cos 18 X"' 
o 


'J n(x")dx" 


( 156 ) 
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vhich, curiously enough, has the same form as the equation given by 
Dungey (loc. cit., page 33) hut is of vectorial character of and includes 
an additional term 3 ( 1 + 5 tan^ x)sec X; the variables and assumptions 
used to arrive at this result differ ratically from those used by Dungey. 

The quantity j(\) = jj_(\)sec3 X verifies, of course, the same 
equation (156). Equation (156) maybe integrated by successive approxi- 
mations . 


Syl vania Electric Products, Inc. 
Waltham, Massachusetts 
May 1 , 1962 
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